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Abstract 

We address the problem of calculating momentum distributions of particles emerg¬ 
ing from the three-body decay of a many-body resonance. We show that these 
distributions are determined by the asymptotics of the coordinate-space complex- 
energy wave-function of the resonance. We use the hyperspherical adiabatic expan¬ 
sion method where all lengths are proportional to the hyperradius. The structures 
of the resonances are related to different decay mechanisms. For direct decay all 
inter-particle distances increase proportional to the hyperradius at intermediate and 
large distances. Sequential three-body decay proceeds via spatially confined quasi¬ 
stationary two-body configurations. Then two particles remain close while the third 
moves away. The wave function may contain mixtures which produce coherence ef¬ 
fects at small distances, but the energy distributions can still be added incoherently. 
Two-neutron halos are discussed in details and illustrated by the 2 + resonance in 
6 He. The dynamic evolution of the decay process is discussed. 


PACS: 21.45,+v, 31.15.Ja, 25.70.Ef 


1 Introduction 


Nuclear resonances and excited states can be very complicated many-body 
structures with a number of different decay modes. The simplest decay, per¬ 
haps beside 7-emission, is breakup into two particles as exemplified by nucleon- 
and a;-emission, and binary fission [1], The deceivingly simple breakup into 
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three particles is much less studied and far from understood. This is in con¬ 
trast to bound three-body cluster structures where a variety of techniques are 
available and able to predict the properties, even of the exotic quantum halo 
states [2], The three-body continuum properties are less established although 
rather well studied over many years [3]. 

Experimental information is obtained by measuring the properties of the par¬ 
ticles in the final state. The experimental techniques are now advanced to 
a level where accurate and kinematically complete measurements are avail¬ 
able on a number of different systems [4-15], and many more are expected 
to follow. Reliable theoretical descriptions are needed to interpret existing 
data, to predict unknown decay results and to help in the design of new in¬ 
teresting experiments. Both the structure of the initial state and the decay 
mechanism are essential and both must therefore be properly described simul¬ 
taneously. Clearly for a genuine many-body state only the intermediate and 
large-distance structure is decisive. The small distance behavior is artificial 
and only serving to provide the proper continuous boundary conditions. 

For two-body decay, like a-emission, the relative potential determines all prop¬ 
erties. In the example of a-emission, the two-body potential can be divided 
into short-, intermediate- and long-distance. The short-distance part is arti¬ 
ficially adjusted to allow the correct resonance energy and the barrier region 
then determines the width. At large distances where the potential has van¬ 
ished the energy of the a-particle is determined by energy conservation. In 
two-body decay the energy distribution then only reflects the width of the ini¬ 
tial resonance. These properties are very different for decays with more than 
two particles in the final state. 

For three-body bound states and resonances “large distance” is less well de¬ 
fined . However, the corresponding structure can efficiently be computed by use 
of the hyperspherical adiabatic expansion method [16,17]. The wave functions 
are in this technique expanded on basis states related to adiabatic potentials 
calculated as functions of the hyperradius p. Then p provides a measure of 
distances for the three-body problem. The wave function is usually dominated 
by the component related to the lowest adiabatic potential. The small-distance 
part (p small) of both wave function and potential are only directly meaningful 
if the particles appearing in the final state form a genuine three-body system. 
Otherwise this part of the effective potential is constructed to produce the 
correct resonance energy and provide an appropriate boundary condition for 
the wave function. At intermediate distances the potential has a barrier which 
is decisive for the width of the resonance. At large distances the resonance 
wave function is characterized by outgoing waves which contain information 
about distributions of relative energies and possibly other quantities like spin 
distributions. 
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For three-body decay the two most obvious decay mechanisms, direct and se¬ 
quential, were recently studied in a schematic model in [18]. Another schematic 
model is also formulated in the limit where only the Coulomb interaction is 
important at intermediate and large distances [19]. The detailed resonance 
structure at small and intermediate distances were investigated in realistic 
models in [20]. The large-distance properties are much more difficult to cal¬ 
culate accurately, because either the correct continuum three-body Coulomb 
wave functions are unknown, or the short-range potentials may produce an 
almost long-range (inverse square) effective potential at large distance. In the 
latter case the origin is precisely as for the Efimov effect [16,21], The corre¬ 
sponding potential is most likely the lowest at large distance but not neces¬ 
sarily also at small distances. Thus, the resonance wave function may change 
structure from small to large distance. The relative energy distribution arises 
as the result appearing at large distance. The numerical computations are 
then sometimes rather tricky. 

The purpose of the present paper is to establish a general method to compute 
relative energy distributions after decay into three particles. The short- and 
intermediate-distance resonance structure from [18,20] is a good starting point 
but we need in addition to calculate the asymptotic behavior in momentum 
space. The asymptotics vary for the different decay mechanisms, and the re¬ 
lated numerical treatment is difficult when all the possibilities simultaneously 
have to be accounted for. 

We assume that formation and decay of the resonances are independent pro¬ 
cesses. The resonance could be formed by beta-decay from a neighbouring 
nucleus, or a window with the relevant energies can be selected in an experi¬ 
mental setup where contributions from other processes also are eliminated. In 
section 2 we develop the theoretical formalism for resonance decay. This was 
previously sketched by use of the saddle point approximation [22], while we 
here shall instead use the Zcldovich regularization of the divergent Fourier in¬ 
tegrals [23]. In section 3 we discuss some of the important features arising from 
calculation of resonance wave functions by use of the hyperspherical adiabatic 
expansion method combined with the complex scaling method. In section 4 we 
illustrate in details with realistic computations for the 2 + -resonance in b He. 
Finally, section 5 contains a brief summary and the conclusions. 


2 Theoretical formulation 


We assume that the system of particles has been generated in a meta-stable 
quantum state (a resonance) that is a generalized eigen-state of the Hamilto¬ 
nian with complex energy. This is a decaying state - it describes a constant 
flux of particles towards infinity. Suppose we have a system of detectors at 
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large distances which measure the momenta of the particles emerging from 
this decaying state. Clearly these detectors will measure the probability dis¬ 
tribution of particle momenta in the meta-stable state, that is the absolute 
square of the momentum space wave-function. 


2.1 Two-body resonances 


The theoretical derivation is most easily understood if we first explain the 
idea for simple resonance decay into a two-body system. We need resonance 
inventions in coordinate and momentum-space and transformations between 
these non-square integrable functions. 


2.1.1 Resonance wave functions 


The momentum space wave-function of a resonant state with the complex 
energy E r = Tj^kf — P 0 — has the form [24] 


tM k ) 


g(k, 0*,) 

1/2 _ 1/2 
lx rv q 


( 1 ) 


where k is the relative momentum and findicates the two angles defining the 
direction of the vector k. We assume that the wave function ipk 0 { k) only has 
the pole at k = k$ and g(k, is then a continuous function of the momentum 
k with no poles. 


The distribution P(k) of the relative momentum k of the two particles in the 
resonant state V’fco is given by the absolute square of the momentum-space 
wave function, i.e. 


P(k) 


l ^ fco ( k )| 2 oc 


\g(k,n k )\' 2 

(P - P 0 ) 2 + ? ’ 


( 2 ) 


where the real observable energy P is P = h 2 k 2 /(2m). The approximation 
that the system is generated in a pure resonant state Vp 0 i s m ost likely only 
valid in the neighborhood of the resonant energy, i.e. P ~ P 0 . Furthermore, 
the function g(k, Tl k ) is smooth and varies by definition much less than the de¬ 
nominator in eq.(l). In any case we shall only consider energies where |P — Pq| 
is less than a few times T. We can then confidently substitute the momentum 
k with the resonant momentum in g(k, Q k ) an d thus arrive at the expression 
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of the famous Breit-Wigner type 


P( k) oc 


|^( fe 0 , | 2 

(■ E-E o y + ^ ’ 


( 3 ) 


where the energy dependence is given by the factor ^(E — E 0 ) 2 + while 

the angular dependence is given by the (absolute square of the) function 
g(ko,Qk)- Thus, the momentum-space wave-function eq.(l) of the resonance 
allows direct calculation of the momentum distributions of the decay frag¬ 
ments through eq.(3). Clearly improvements are possible by use of different 
approximations of g(k, Q* ; ). 


Instead of momentum-space, the wave-function of the resonance may be avail¬ 
able in coordinate-space where the large-distance asymptotic form is given 
by 




D +ik 0 r 


'/(fir) , 


( 4 ) 


where fl r denotes the two angles defining the direction of the relative coordi¬ 
nate r. The structure is an outgoing spherical wave potentially modified by an 
angular dependence contained in /(h2 r ). Generally the resonance wave func¬ 
tion can be written as a partial-wave expansion in the spherical harmonics 
Y tm , i.e. 

fM r ) = J2 C lmXl(r)Ylm(^r) , (5) 

lm 


where C) m are constants depending on angular momentum and projection 
quantum numbers / and m. The radial functions Xi( r ) are those solutions of 
the radial Schrodinger equation that asymptotically approach the outgoing 
spherical wave in eq.(4), i.e. 

e +ik 0 r 

Xl(r) ™ - , /(fir) = E ClmYlm^r) ■ (6) 

r , 

lm 


This defines the asymptotic behavior of the decaying resonance wave function 
which in turn determines the energy distribution in the observable final state. 
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2.1.2 Transformation from coordinate- to momentum-space 


The coordinate-, tfk 0 ( r )> and momentum-space, ifk 0 (k), wave-functions are 
connected via a Fourier transform 

f’ko ( k ) = f e~ tkr if ko (r)d 3 r. (7) 

Expansion of the plane-wave in terms of spherical harmonics 

e ikr = £ 4 tt i'jiWYUQjYtKQk) (8) 


reduces the Fourier integral in eq.(7) to a one-dimensional radial integral 

OO 

^ fco (k)=47T YhCim{-i) l Yim(£lk) f r 2 drxi{r)ji{kr) . (9) 

lm n 


Because of the asymptotics in eq.(6) the radial integral is seen to diverge. The 
large-distance behavior is responsible for the divergence. The physics content, 
expressed by a finite value, then has to be extracted by a suitable regulariza¬ 
tion. We use the prescription proposed by Zcldovich [23], i.e. multiplication of 
the integrand by a gaussian where the range after integration is increased to 
infinity. For an exponential this gives 

OO OO _ _ 

[ e iqr dr -> lirn / e iqr ~ a2r2 dr = lirn erfc(-—) = - , (10) 

1 <*->o J a->o 2a 2a q 

0 0 


where q can be complex and erfc is the error function of complex argument. 
In the present context the radial integral in eq.(9) is first rewritten by adding 
and subtracting the asymptotic expression of the diverging integrand. The 
difference between the true and the asymptotic expression then remains finite 
even without multiplication by the gaussian function. Only the asymptotic 
expression then diverges when the gaussian smoothly converges to an overall 
factor of one. 

The physics content is extracted by dividing with a similarly diverging nor¬ 
malization integral of the square of the wave function Xi( r )- Also this integral, 
now in the denominator, is rewritten by addition and subtraction of its asymp¬ 
totic expression. Again only the asymptotic expression diverges. The Zcldovich 
prescription now leaves the ratio of these two diverging integrals of the asymp¬ 
totic expressions. However, this ratio does not diverge but converge towards 
the physically meaningful result. Apart from a normalization we therefore have 
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to regularize only the asymptotic expression obtained by use of eq.(6) and the 
asymptotic approximation of ji(kr), i.e. 


V’fco (k) = 47 T^2Ci m (-i) l Y lm (tt k ) 

lm 



e +ikor 


r 


sin (hr — y) 
kr 


( 11 ) 


The radial integral is then by use of eq.(10) evaluated to be 


OO 



0 



2 


1 

k — k 0 


i l k + k 0 - (-1 )'(k - ko) _ k 0 

o’ U'2 _ U2 — 1 U2 _ U2 

Zj IXj /VQ IV Iv Q 


(- 1 ) 1 ' 
k + ko 

■i k 

_ 

U 2 _ 1,2 

rv A.Q 


( 12 ) 


for even or odd l, respectively. The summation in eq.(ll) is proportional to the 
angular amplitude / from eq.(6), but now as a function of the momentum k. In 
any case we assumed earlier that k ~ ko in the smooth functions. We therefore 
arrive at the final expression for the Fourier-transform of the resonance wave- 
function, i.e. 

Ajr _ Ajr 

<Mk) = Y 2 72 E C, m Y lm (Sl k ) = — p/(Sl t ) , (13) 

K K 0 lm K K 0 


Thus the function g from eq.(l) is then related to / from eq.(4) by 


g(k 0 ,tt k ) = 4vr/(f] fc ) . 


(14) 


The convenient fact that only the asymptotic limit of the resonance wave 
function enters after the regularization procedure is perhaps more surprising in 
mathematics than in physics where the observable energy distributions always 
are obtained from the properties at large distances. 

The observable distribution in momentum-space is determined by the angular 
wave function in coordinate-space evaluated for angles describing the direction 
of the momentum. This peculiar fact can intuitively be understood by the ge¬ 
ometry of particles moving towards the detectors at infinitely large distances. 
Coordinates and momenta then must point in the same direction. A mathe¬ 
matical formulation is available from ionization cross sections calculated for 
atomic physics processes [25]. 
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2.2 Three-body resonances 


The generalization to three particles first requires a convenient set of coordi¬ 
nates. We choose the scaled Jacobi coordinates [16] 


x = 


m 2 m 3 


m(m 2 + m 3 ) 


( r 2 - r 3 ), 


y = 


N 


mi(m 2 + m 3 ) 
m(mi + m 2+ m 3 ) 



m 2 r 2 + m 3 r 3 \ 

m 2 + m 3 J ’ 


(15) 


where m is a mass scale, and r., and rri t are the coordinate and mass of the 
particle number i. The hyper-spherical coordinates are then the hyper-radius 
p, the hyper-angle a, and the directional angles Q x and VL y of the vectors x 
and y defined by 

p = \Jx 2 + y 2 , a — axctan(x/y) , Q p = {a, Q x , . (16) 


The Jacobi coordinates depend on the sequence chosen for the particles, and 
the three different pairs of x and y could be labeled to distinguish. We omit 
these labels when the meaning is clear. 

The corresponding momentum-space variables are 

k = \jv 2 + q 2 , « K = arctan(g/p) , f \ , (17) 


where p and q are the conjugate momenta related to the coordinates x and 

y- 

2.2.1 No bound two-body subsystems 

The generalization of the two-body spherical harmonics are the so-called hyper- 
spherical harmonics [16] 


yic(^p) = N% x,ly ^ sin*® a cos ly aPn X+ 2 ’ ly+ 2 \cos2a) 

XYl*m x (ttx)Yl y m y (n y ) (18) 

where /C = {Kl x m x l y m y }, K = 2n + l x + l y , and (l Xl rn x , l y , rn y ) are the angular 
quantum numbers related to coordinates x and y, and N^ 0,ly ' ) is a normaliza¬ 
tion factor. These functions are the eigen-functions of the angular part A 2 of 



the three-body kinetic energy operator T 


T 


2m 


iyl + 


n 2 \_d 2 _ 5 d A 2 

2m dp 2 p dp p 2 


(19) 


with the eigenvalues K(K + 4), i.e. 


A 2 y lc (n p ) = K(K + 4)y lc (n p ), 


( 20 ) 


where K is a non-negative integer. Without Coulomb and without bound two- 
body subsystems the three-body resonance wave-function T K0 (p, Q p ) with the 
complex energy E r = H 2 k,q/( 2 m) = E 0 — ?T 0 /2 can be expanded in terms of 
the hyper-spherical harmonics 

{Pi Qp) = 5Z C/cXic(p)yic(ttp)i (21) 

K 


where the hyper-radial functions Xic(p) have the usual resonance asymptotic 
behavior of an out-going hyper-spherical wave 


xM - 


g+lKOP 

"yv- 


( 22 ) 


The three-body wave-function asymptotically has the form of the hyper-spherical 
wave with an angular amplitude F(Q P ) determined by the expansion coeffi¬ 
cients Cfc, i.e. 


pTtrvup_ p-ror^u/j 

*«.(/>, n„) ^ -ITT E = -Jji (23) 

P K P 

The momentum-space wave-function is the Fourier transform 

t ko («, n K ) = I e - ipx - iqy T K0 (p, n p )p 5 d P dn p . (24) 


The three-body plane-wave can be expanded in hyper-spherical harmonics as 


\ K P) K. 


(25) 
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Due to orthogonality of the hyper-spherical harmonics the angular part of the 
integral in eq. (24) is trivial and we are only left with the hyper-radial integral 


\I/ K q(^, 0 K ) — ^ } y ^(—i) K Cjcyic(Q K ) [p 3 dpxic(p)JK+2(^p)- (26) 

h\i i^ J 


Precisely as in the two-body case, in the vicinity of the resonance the inte¬ 
grand can be made ready for regularization by substitution of eq.(22) and the 
asymptotic form 


Jk+2(kp) 


p—XX) 


2 irK 

-sm (up - — . 

7 xnp 2 


(27) 


This results in a diverging integral similar to that of the two-body case 

**,(«, OO = E(-*) S 'C K y K (S2«) jdpe+^sinUp- ^).(28) 

' /c 


Using the Zcldovich regularization leads to 


Kq V tt ^ z a. — h , 0 

27/2^5/2 X 

= 372 77 (29) 

K 0 K — Kq 

that is, in the vicinity of the resonance, the angular wave function in momentum- 
space is proportional to that in coordinate-space but evaluated for the mo¬ 
mentum variables. The energy distribution is determined by the Breit-Wigner 
factor where the width is obtained from the three-body resonance. The func¬ 
tion F(Q k ) now contains information about the non-trivia! energy distribution 
between the three particles. This is in contrast to the two-body decay where 
all the energy is in the only existing relative degree of freedom. 


2.2.2 One bound two-body subsystem 

Sometimes a bound two-body subsystem can be emitted from a three-body 
resonance. Such a final state configuration can not be described by hyper- 
spherical harmonics. However if this is the only open channel, the description 
of such a decay reduces to the two-body case. Indeed the asymptotics of the 
resonance wave-function is then 

P my 

T K0 (p, Q p ) ^ fc3( x ) /(fy), (30) 
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where 023(x) describes a bound system of particles 2 and 3 with binding 
energy B 2 3 , qo = \j2m{E r — B 23 )/h 2 , f(0 y ) is the angular amplitude and E r 
is the complex three-body energy. If both three-body and two-body decays are 
possible the wave-function contains asymptotics of both two- and three-body 
types, 


(P) ^p 


p—XX) 


A 


e +iK 0 p 




P 


,5/2 


F(fi p ) + 50 23 (x)— f(n y ), 


(31) 


where A and B are the asymptotic coefficients determining the relative weights 
of the two decay channels. Both F and / are dimensionless and the dimension 
(length to —3/2) of 0 2 3 compensate for the one length in the denominator of 
the last term. 

The Fourier transform and the corresponding regularization then give the 
momentum-space wave function, i.e. 

97/2 5/2 a 4 

^«o( K > ^k) —--2- 2 F (^ k ) + B ( p 23 ( p )— - 2 fy (^ q ) j ( 32 ) 

Kn ~ K ~ K 0 Q ~ % 


where q 2 = k 2 — 2mB 23 /h 2 , 023 (p) is the momentum-space wave function of 
the two-body bound state 0 23 . These two channels correspond, respectively, 
to two close-lying particles in a bound state far away from the third one, and 
three particles all far away from each other. Thus, in the of large distances 
they do not interfere, and the resulting momentum distribution is simply a 
weighted sum of the corresponding distributions. The relative contributions 
of the two channels are given by \A\ 2 and \B\ 2 correspondingly. 

Generalization to describe decays into more than one two-body bound state 
is formally straightforward, i.e. the corresponding non-interfering asymptotic 
terms should simply be added. This holds for more than one bound state in 
the same two-body system as well as for bound states in different two-body 
systems. 


2.2.3 One resonant two-body subsystem 

Instead of a bound state the decay via a two-body resonance is often con¬ 
sidered in interpretation and analysis of experiments. Clearly the narrower 
the resonance the more similarity to the case of two-body bound states. In 
any case the hyper-spherical expansion must eventually converge, although 
the convergence can be too slow for a reliable extraction of the asymptotic 
coefficients in eq.(23) from a numerical solution of the three-body problem. 
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However, in this case the (slowly convergent) two-body resonance configu¬ 
ration can then be explicitly included into the asymptotics while only the 
remaining (hopefully fast convergent) part is expanded, i.e. 


^{p,n p )^A 


e +iKop 


D ipoX 




P 


5/2 


F{n p ) + B — f x {n x ) — fy(Qy), 


X 


y 


(33) 


where h 2 K^/{2rn) = Eq—iTq/2 is the complex three-body energy, h 2 pl/(2m) = 
E 23 ~ ^ 23/2 is the (complex) energy of the two-body resonance and the 
remaining part is described by the complex momentum q g = — p 2 . The 

precise definition of arises from a constraint to be seen below. 


The corresponding momentum-space wave-function is again given by the reg¬ 
ularized Fourier transform, i.e. 


T K 0 (k,H k ) = -A 


27 / 2 ^ 5/2 


K, 


3/2 


+ B- 


47T 


K 


Kn 


pz 


Po 




47T 


V 9 

qi 


% 


m,). ( 34 ) 


The momentum distribution is given by the absolute square of this momentum- 
space wave-function. In the center of mass system we can directly find the 
distribution of particle 3 arising from the sequential emission via the two- 
body resonance. 

Absolute square of the last term in eq. (34) and use of the energy conservation 
k 2 = q 2 + p 2 (or E = E\ + E 23 ) immediately gives the energy distribution for 
particle 1 


P(E!) oc 


dEa |(S 23 - -E#) 2 + rls/4] pia + El - Bo) 2 + rg/4] 

1 

(X - 

(El - (E 0 E ^) 2 + (To + r 23 ) 2 /4 ’ 


(35) 


which states that the most probable energy of particle 1 is E\ = Eq — E^ 
and the width of the distribution is the sum of the two and the three-body 
widths. Precisely this Breit-Wigner distribution only arises when the \q 2 — q ^ 2 
in eq.(34) is proportional to ( E — Eq ) 2 + Tq/ 2, i.e. given by the probability 
distribution in the initial three-body state, which also is of Breit-Wigner form. 
Thus the definition of must involve p 2 and not Pq. 


The same integration could of course be performed on the energy of particle 
1, but this would only give the two-body Breit-Wigner distribution for A 23 
whereas the measurements provide individual energies, E 2 and E :i . for particles 
2 and 3 in the center of mass system. To get the distributions of E 2 and E :i 
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involve trivial but tedious kinematical transformations where also energies and 
directions of particle 1 are required. 


2.2.4 Alternative real-coordinate procedure 

The relative energy distributions can be obtained from the angular resonance 
wave function calculated without complex scaling. We assume that the system 
of particles is produced in an initial state for example by a beta-decay process. 
We can then imagine the subsequent decay as the time evolution of the initial 
non-stationary state. This can be formulated as a time dependent coupled 
channels problem. It can also be viewed intuitively as a particle described 
by time dependent coordinates determined by classical equations of motion. 
This should be done with the appropriate initial amplitudes for all parts of 
the initial wave function. The hyperradius must vary from being very small 
to infinitely large. We increase p until all particles are outside the interaction 
ranges of all other particles. From then on all distances scale as p and all other 
coordinates remain unchanged with time until p = 00 . 

Energy conservation is maintained at large distances by converting the po¬ 
tential energy into kinetic energy in the scaling degree of freedom, i.e. by 
increasing the velocity p of the p-coordinate. The wave function then evolves 
with all angular degrees of freedom frozen eventually reaching the detectors 
placed far away. The absolute square of the angular wave function as function 
of cos 2 a then provide the energy distributions simply because the kinetic en¬ 
ergy of particle 1 is given by the velocity of y , i.e. y = p cos a. Then the energy 
distribution as function of the kinetic energy, proportional to y 2 oc cos 2 a, is 
the probability coordinate-space distribution as function of cos 2 a, apart from 
the phase space conversion from a to energy, i.e. division by a factor propor¬ 
tional to dE/da oc sin(2a). 

This procedure is tempting since we only need to increase the maximum value 
of p in all the numerical implementations and plot the wave function at that 
large distance. For this to be accurate the asymptotics has to be well described 
by the hyperspherical expansion or the basis has to be very large and able to 
describe the necessary large p-behavior. However, when the basis functions 
have asymptotics different from one of the intermediate structures the size of 
the basis needed to reach convergenge can easily become huge, making the 
procedure impractical. This is not necessarily easy to see in the numerical 
results where an increase of basis size usually is rather expensive while the 
convergence could be extremely slow. The procedure is probably only directly 
useful for direct decay or for sequential decay via broad resonances. Other¬ 
wise the different intermediate structure should be computed somehow and 
extrapolations to large distances applied to each component individually. 
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An example to illustrate the present alternative formulation is available in 
the schematic model discussed in details in [18] where the widths but not the 
energy distribution were computed. Assume that only the Coulomb potential 
is important and the most probable path (ridge in the wave function) from 
small to large distances can be described by scaling the hyperradius. The 
corresponding optimum path is defined by minimizing the WKB-tunneling 
expression as function of different relative scaling parameters Sik = r^/p, 
where r lk is the distance between particles i and k. The path is given by 
sf k rriiZj = s^ k rrijZi (see also [19]), where Zje is the charge of particle i. We 
then arrive at the most probable value for the energy division, i.e. 


Ek 



™> k Zl i/ 3 \ 1 
rnjZy ) ’ 


(36) 


where E tota i is the total energy distributed among all the three particles. This 
expression is simple but not very accurate, ft also only provides an estimate of 
the peak value. To compute the distribution other paths must also be consid¬ 
ered. This is possible but we shall leave this for a later discussion in connection 
with a detailed treatment of the Coulomb interaction. 


3 Resonance wave functions 


The resonance wave function contains all information including that of the rel¬ 
ative energy distribution after the decay. The calculations must then first pro¬ 
vide the corresponding three-body resonance states. Second the large-distance 
behavior must be accurately extracted. Due to the different structures, these 
steps are not trivially completed by use of only one method. We briefly de¬ 
scribe first the main ingredients in our computations and the features of the 
wave function. Second we explain how the asymptotic behavior is obtained in 
practice. 


3.1 Method 


We use the hyperspherical adiabatic expansion method combined with com¬ 
plex scaling to obtain resonance wave functions. The coordinates are defined 
in section 2 along with our basis functions in angular space, i.e. the hyper¬ 
spherical harmonics in each of the three Jacobi systems. We solve the com¬ 
plex scaled Faddeev equations as function of hyperradius [16,26]. The complex 
scaled coupled set of radial equations are subsequently solved with the appro¬ 
priate boundary conditions, i.e. exponentially vanishing with increasing p for 
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both bound states and resonances. Thus here we assume that we do not need 
to treat the Coulomb interaction explicitly at asymptotic large distances. As 
shown in [26], the results obtained with this method agree well with some of 
the most common procedures, as for instance the complex energy method. 

The energies are usually accurately determined in this method. The same 
applies to the wave functions at small and intermediate distances where the 
exponential fall-off still is not too restrictive. However, we need the informa¬ 
tion at distances where the asymptotic limit is reached, i.e. possibly at very 
large p where the complex scaled wave functions are very small. Furthermore, 
more than one geometric structure can be important at the same time, e.g. 
two different spatially confined two-body configurations with the (different) 
third particle far away. This happens frequently with two identical particles 
like neutrons and protons as constituent particles, because the nucleon-core 
interaction must be sufficiently attractive to produce a bound or resonating 
three-body system, and this implies that such two-body configurations are 
favored. To account simultaneously for different two-body substructures, it is 
essential to use three components as in the Faddeev decomposition adopted 
by us. The same efficiency can be achieved in a variational approach by allow¬ 
ing Faddeev-like components in the trial wave function [27]. It is much more 
difficult, if not impossible, to reach convergence with only one component as 
in the hyperharmonic expansion method [28]. 

Even with three Faddeev components a large basis has to be employed. To 
describe substructures inside one of the two-body potentials (range R e ff) for 
large p we need values of the hyperpsherical quantum number K up to a 
few times P\frnj (R e ff^p) where p is the reduced mass of the two particles. 
This is because K/2 is the number of nodes in the basis, and details can only 
be described if a few nodes can be placed inside the structure in question. 
Thus, for nuclear systems, where R e ff ~ 4 fm, we need K max of at least 
50 to describe such structures for p ~ 100 fm. Employing complex scaling 
transforms resonances into states obeying the numerically easier bound state 
boundary conditions. The required basis size is larger, essentially because the 
exponential fall-off moves to larger distances with increasing scaling angle. 
These estimates provide necessary conditions for a reasonable description of 
two-body substructures which in turn are necessary to describe the sequential 
decay mechanism. 

The requirement of large K to describe substructures must be reconciled with 
the fact that an increase of p towards infinity results in convergence of the an¬ 
gular eigenvalues to the hyperharmonic spectrum for free particles. We show 
an example in Fig. 1 where the imaginary parts are omitted as they both 
oscillate around zero and approach zero at large p. The real parts of the an¬ 
gular eigenvalues approach K(K + 4) as p increases while the corresponding 
potentials all approach zero faster than p ~ 2 . The attractive pockets in the 
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Fig. 1. The real parts of the lowest 8 angular eigenvalues (left) and corresponding 
adiabatic potentials (right) as functions of p for the 2 + states in 6 He ( 4 He + n + 
n). The scaling angle is 6 = 0.10 rads. 

eigenvalues at short distance disappear in the potentials except for the two 
lowest where the negative values remain as a prominent feature. The approach 
to the asymptotic values is very fast except for the levels where s-waves con¬ 
tribute significantly. The low energies favor these levels at large distances. 

The related adiabatic wave functions approach the hyperspherical harmonics. 
The reason is that the regions in space, where the short-range interactions are 
significant, are shrinking in size with increasing p relative to the total space 
available. The radial extension of these regions, responsible for two-body corre¬ 
lations, decrease with p as 1/p. The interactions are non-vanishing in smaller 
and smaller regions. Consequently they become less and less important for 
both energies and wave functions. Thus, the basis size has to increase with 
p in order to allow a description of the two-body substructures or equiva¬ 
lently of sequential decay, but the lowest adiabatic potentials approach the 
free solutions. The basis size in practice always has to remain finite and the 
substructures eventually become impossible to describe in this way. 

The interactions have to be chosen to reproduce the pairwise low-energy scat¬ 
tering properties of the three particles appearing in the final state. Clearly we 
must accurately include all the partial waves necessary to describe the quan¬ 
tum numbers of the decaying resonance. However, even with a sufficiently 
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large basis the three-body system does not necessarily appear with the cor¬ 
rect energy and width. In fact, there may not even be an attractive region at 
small distances as required to produce a resonance of finite width. This could 
occur when we are dealing with a many-body resonance without traces of any 
three-body cluster structure. Nevertheless, a meaningful computation can be 
carried out of the energy distributions emerging after a three-body decay. 

The philosophy is the same as for ct-emission where the inner part of the effec¬ 
tive potential is replaced by an attractive square well with a depth adjusted 
to reproduce the resonance energy. The resulting barrier is then used to derive 
the width, usually in the WKB approximation. We generalize this concept to 
the adiabatic potentials, i.e. we add a three-body potential of short range in 
the hyperradius. It is intended to describe interactions beyond the two-body 
level such that the three-body system has a resonance at the desired energy. 
By doing this we have substituted the possibly complicated many-body struc¬ 
ture at small distance with the three-body cluster structure resulting from an 
effective potential, which in turn also provides the correct boundary condi¬ 
tions for a three-body decaying resonance. This principle was introduced for 
fine-tuning in the first calculation with the correct boundary conditions of 
the three-« decay of the second 0 + -state in 12 C [29]. It has later become the 
standard procedure to adjust three-body energies without significant changes 
of the underlying substructure [16]. 


3.2 Important features 


The radial solution is often strongly dominated by one or two of the lowest 
adiabatic components at small distance where the relative probability is large. 
This is because all three short-range two-body interactions contribute simulta¬ 
neously and the result is the energetically most favored three-body resonance 
structure consistent with the boundary conditions. As p increases, at least one 
particle has to move away from the other leaving at most one non-vanishing 
two-body interaction. Coherent contributions from different of these config¬ 
urations are possible and sometimes even favored. At large distance, where 
the energy distribution is determined, several more adiabatic potentials are 
often needed. The couplings due to the Coulomb interaction would normally 
increase the necessary number of potentials. 

It is established [16] that the adiabatic potentials of lowest energy at large 
p are related to configurations with relative s-waves between the two closest 
particles. This is the basis for the Efimov effect [30]. If these large-distance con¬ 
figurations differ from the resonance structure at small p, possibly with higher 
partial waves, the lowest angular wave function must change its structure ac¬ 
cordingly as p increases. In [20] we showed the structure for 6 He(2 + ) for the 
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Fig. 2. The fraction of different components in the fifth adiabatic potential for 
6 = 0.10 rads as function of p for 6 He(2 + ). The angular eigenvalue corresponds to 
K = 4 at large p, see Fig. 1. The angular momenta are specified by £ x , j x , £ yi j y , 
and L. Left: x refers to the two-neutron system and y to its center of mass motion 
relative to the a-particle. Right: x refers to the neutron-ct system and y to its center 
of mass motion relative to the other neutron. We give the (,x, y) components on the 
figure as £j. 

lowest eigenvalue which approach the K = 2 value at large p. In Fig. 2 we show 
the results for the similar eigenvalue approaching the K = 4 level for large p. 
The pronounced and rapidly changing structure is qualitatively similar to the 
lower-lying K = 2 level, i.e. dominated by p 3 / 2 — p 3 / 2 neutron-core structure 
at small p and by s-waves between the two neutrons at large p. Essentially all 
other allowed components contribute with equally small amounts. 

These rather dramatic changes imply that it is crucial to include all adia¬ 
batic potentials with significant couplings to those dominating the structure 
at small hyperradii. This is simply because the couplings are responsible for 
changing the radial weights of the different adiabatic components as func¬ 
tion of p, e.g. no couplings imply the same occupation independent of p. On 
the other hand, each of the angular wave functions related to the adiabatic 
potentials are themselves functions of p, sometimes rapidly varying as seen 
in Fig. 2. In principle the non-diagonal couplings could be vanishingly small 
and all change of structure would be described by the lowest adiabatic wave 
function. However, this is rather unlikely because the couplings are defined 
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Fig. 3. The coupling potentials between the four dominating adiabatic levels for 
9 = 0.10 rads shown as functions of p for 6 He(2 + ). The first and the fourth levels 
have similar quantum numbers but approach the K = 2 and 4 levels, respectively. 
To show the first (P) and second (Q) order coupling potentials in the same units 
(fm -1 ) we multiply Q by p. (The energy unit is restored in the coupling potentials 
by including the omitted factor, i.e. h 2 Q/(2m), h 2 P/ (2m)d/dp). 

as matrix elements of first and second radial derivatives of the angular wave 
functions. Thus, radial couplings between rapidly changing angular structures 
are inevitable. 

In Fig. 3 we show how the strongly varying coupling potentials can be re¬ 
lated to the changing angular structure seen in Fig. 2. The peaks are most 
pronounced when a crossing is avoided and the two levels switch characteris¬ 
tics [16]. The rather confusing coupling picture is crucial for the asymptotic 
behavior of the wave function at large distance where the energy distribution 
is determined. The second order terms are substantially larger than the first 
order couplings but all vanish at large p. Thus the numerical computations in 
the present case must extend at least beyond p ~ 40 fm where the couplings 
have reached very small values. This is about two times the largest scattering 
length which defines the distance of convergence towards asymptotic values. 

We have now established two important but competing effects, which deter¬ 
mine the three-body resonance structure from small to large values of p. In the 
extreme, the structure can either remain unchanged by climbing correspond- 
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ingly up on the adiabatic potentials, or the structure can change to follow that 
of the lowest-lying adiabatic wave function. A compromise between following 
the energetically most favored configuration and the resistance to a change of 
structure therefore must be reached. The combination of these effects deter¬ 
mine the relative population of the different components in the radial solution, 
which in turn determines the observable energy distribution. The couplings are 
more important here than for widths, energies and small-distance wave func¬ 
tions. They must be accurately computed to provide the energy distribution. 

The structure of the resonance wave function at large p could remain un¬ 
changed and only exhibit a simple scaling behavior proportional to p. This is 
typical of direct decay. The wave function could also have large probability 
for finding two close-lying particles where the hyperradius mainly changes by 
moving the third particle as p increases. This is typical of sequential decay via 
more or less stable two-body configuration, e.g. sequential decay through two- 
body resonances. The intermediate configuration does not necessarily need 
a confining barrier, but could be provided by low-lying two-body virtual s- 
states [21]. Mixtures of all types can occur giving rise to the description of 
decay properties as fractions proceeding via individual two-body configura¬ 
tions. All these structures can be computed by use of our method, although 
convergence for the Coulomb interaction is more difficult. 

The best strategy to get reliable results is not obvious, because the brute force 
method of increasing basis size and hyperradius until convergence is reached 
may be beyond any reasonable computer effort. The indecision is related to 
the requirement of an increasing basis with increasing p, which means that a 
smaller p and a smaller basis could provide a better description with much less 
effort. In other words a convergence may be reached in a region of p -values 
for a moderate basis size. This convergence would be destroyed as p is allowed 
to increase because the basis size cannot follow. The convergence can possibly 
be reached faster by extrapolation of the observable distribution by use of a 
known or anticipated dependence of p and basis size [31]. Different parts of the 
wave function may extrapolate differently. The most efficient choice depends 
on the (mixtures of) decay mechanisms which therefore has to be determined 
first. Therefore the first step is to compute the structure of the resonance wave 
functions as discussed in [20]. 


4 Realistic numerical illustration: 6 He(2 + ) 


Nuclear three-body decay without complications of the Coulomb interaction 
must involve emission of two neutrons. The decaying states do not have to 
be three-body structures although such two-neutron halos are available and 
rather well studied. The most obvious case is the established 2 + resonance in 
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6 He which is formed by the same neutron-core components as in the ground 
state. 

Without Coulomb interactions the computations should quickly lead to the de¬ 
sired energy distributions. However, even short-range interactions can present 
difficulties as highlighted by the intricate description needed for the Efimov 
effect [16]. Both the a-neutron and the neutron-neutron interactions are pre¬ 
viously employed in ground state computations [16]. 

We follow the procedure outlined in the preceding sections. Different prescrip¬ 
tions are possible to implement the Pauli principle [32], all of them provid¬ 
ing indistiguishable angular wave functions at large distances. We adjust the 
three-body potential to give the correct resonance energy. This only requires 
marginal fine-tuning. In total 1132 hyper-spherical harmonics are used in the 
expansion (21), and the maximum value of K is 200 for the most relevant par¬ 
tial wave components and never smaller than 40. The resonance wave function 
is then available as function of the hyperspherical coordinates. A complex scal¬ 
ing angle of 0.10 rads is enough to produce an exponentially vanishing with 
increasing p wave function for the resonance. A different scaling angle, where 
the numerical calculations have converged, produce the same results. 


4-1 Resonance structure 


We already showed the angular eigenvalues and the adiabatic potentials in 
Fig. 1. The probability distribution arising from only the lowest potential was 
shown in [20]. The structure changes from peaks at small p corresponding to 
a-neutron p 3 / 2 -structure to a probability with one broad peak corresponding 
to comparable distances between all three particles. This reflects the change 
of structure of this angular wave function from small to large p as seen in 
details in Fig. 2. Eventually the lowest hyper harmonic function with K = 2 
is approached. This indicates in itself a direct decay mechanism. However, in 
this case the lowest potential provides rather misleading results. 

The rapidly changing structure seen in Fig. 2 at around p « 20 fin could 
easily lead to occupation of higher-lying levels. These occupation probabilities 
are functions of p and simply found as squares of the radial wave function 
obtained by solving the coupled set of radial equations. The results are shown 
in Fig. 4 for the lowest adiabatic components. At small p the lowest potential is 
totally dominating but as p increases the lowest three components contribute 
with comparable amplitudes. All radial wave functions vanish by oscillating 
around zero with decreasing amplitudes. The relative sizes are more clearly 
seen in the right hand side of Fig. 4. After the transition around 20 fm the 
individually very small radial amplitudes stabilize on relatively constant finite 
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Fig. 4. The radial wave functions (left) and the absolulte values and real parts of 
their relative sizes (right) corresponding to the four dominating adiabatic potentials 
for 6 = 0.10 rads as functions p for the 6 He(2 + ) resonance. 



Fig. 5. The probability distribution for 6 He(2 + ) including the lowest 8 adiabatic 
potentials as function of hyperradius p and hyperangle a related to the distance by 
rik oc ps'ma. i.e. the distance between either the one neutron and core r nc (left) or 
the two neutrons r nn (right). 

ratios. The square of these give the relative weights, i.e. reduced compared to 
the first component by about 0.6, 0.25, 0.01 for the second, third and fourth 
potential, respectively. The transition to stable ratios is consistent with the 
disappearance of the coupling terms shown in Fig. 3. 
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Fig. 6. The energy distributions of neutrons (right) and a-particles (left) after de¬ 
cay of 6 He(2 + ) for 6 = 0.10 rads. The three-dimensional plot show the depen¬ 
dence on p with inclusion of 8 adiabatic wave functions. The maximum energies 
are ( m a + m n )/(m a + 2m n )Eo and 2 m n /(m a + 2m n )Eo for the neutron and the 
a-particle, respectively. Here Eq is the energy of the decaying resonance. 

The total probability distribution in Fig. 5 are quite different from that of 
the lowest eigenvalue. At large p the probability now peaks at a smaller dis¬ 
tance between the two neutrons and correspondingly the a-neutron distance 
is increased. Still fairly broad distributions remain. The decay mechanism in¬ 
dicated by this structure is now instead of direct rather a mixture between the 
preferred sequential decay via a neutron-neutron intermediate configuration 
and a smaller direct component. 


4-2 Energy distributions 


Reliable computation of the energy distribution requires numerically con¬ 
verged results in an appropriate region of p-values. The energy distributions 
are shown in Fig. 6 as functions of p for a sufficiently large number of adi¬ 
abatic potentials. The resemblance with the probability distribution is not 
surprising since only the volume element has been changed. The observable 
distribution is the cut for constant, and sufficiently large, p where convergence 
has been reached as function of basis size. The neutron energy distribution 
has two peaks for small p corresponding to the geometric configurations of one 
neutron close to the a-particle and the other neutron further away. This is re¬ 
flected in the peak in the a-spectrum at intermediate energies corresponding 
to the same geometric configurations. 

The structure changes with p into a broad peak at intermediate energies for 
the neutron spectrum, and one peak very close to the maximum energy for 
the a-spectrum. This is the fingerprint of sequential decay via emission of the 
a-particle followed by decay of an intermediate two-neutron structure. This is 
easily visualized as the two-body decay process where the a-particlc receives 
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Fig. 7. The energy distribution of the a-particle after decay of the 2 + -resonance in 
6 He. The scaling angle is 0 = 0.10 rads and p = 75 frn where convergence is reached. 
The points are extracted from the measurements in [4], Contributions from the 4 
dominating adiabatic potentials are shown individually. 

maximum energy when the two neutrons move together in the opposite di¬ 
rection. In the subsequent decay each neutron then must share the remaining 
energy which leads to an intermediate energy between zero and the maximum 
value. 

This inferred decay mechanism is perhaps counter-intuitive because stable 
intermediate configuration of two neutrons do not exist neither as bound states 
nor as resonances. It would be much more acceptable with the a-neutron P 3 / 2 - 
resonance as the intermediate structure. However, one characteristic feature of 
the neutron-neutron interaction is the low-lying virtual s-state which simply 
means that there is a substantial a-wave attraction. Apparently this is decisive 
for the decay process where the two neutrons end up by moving essentially 
in the same direction, and then necessarily guided by the attraction. The 
interesting point is maybe that this is not the way they started out at small 
distance in the spatially confined part of the wave function. This change of 
structure with hyperradius is in a sense reflecting the dynamic character of the 
decay process. At distances larger than the scattering length the short-range 
interactions are negligibly small, the wave function changes are completed and 
the asymptotics are established. 

The microscopic structure of the energy distributions can be studied by di¬ 
viding into contributions from the individual adiabatic potentials as seen in 
Fig. 7 for the emitted a-particle. The total distribution remains essentially 
unchanged if more than the four dominating potentials are included. Each 
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Fig. 8. The same as Fig. 7 for the neutrons emerging after decay of the 2 + -resonance 
in 6 He. 

contribution has its own characteristic feature. The first has a peak close to 
the maximum energy, i.e. resembling a-emission from a neutron-neutron 0 + - 
state. The second has a peak at intermediate energy, i.e. resembling sequential 
decay by the a-neutron resonance. The third has a peak at small energy, i.e. 
resembling a-emission from an excited neutron-neutron 2 + -state. In addition 
the fourth potential also gives a small contribution with maxima at inter¬ 
mediate and maximum energy. The size of about 1% cannot be seen in the 
total distribution. However, this eigenvalue has the same angular momentum 
quantum numbers as the first level. Therefore the non-diagonal interference 
term would be about 10% of the total contribution. It turns out that the in¬ 
terference essentially is destructive and responsible for the almost flat region 
at intermediate energies. 

The decay mechanism is then not simple although understandable in terms 
of our formulation. The main contribution is decay via the virtual s-state 
and the second is from direct decay. The third mechanism is produced by 
the coupling to the higher-lying state taking place at relatively small p. This 
populates the level eventually approaching the K = 4 hyperspherical level at 
large p. The interference with the dominating contribution then leads to the 
total distribution. The division into these different contributions of direct and 
sequential is to some extent artificial but perhaps useful in connection with 
the experimental analysis. 

The mixture of all these contributions leads to the total distribution which has 
the right features but without precise reproduction of the high-energy peak, 
see Fig. 7. Experimental resolution would not improve very much because ei- 
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ther the peak gets broader and lower, or it gets higher and narrower. The 
discrepancies can originate from the presence of the target and the reaction 
mechanism itself as well as from contributions to the experimental points from 
other than resonance decays. The experiment selects the window of energies 
around the 2 + resonance position in the reaction 'Li( 2 H, 3 He) 6 He*. This nec¬ 
essarily includes some background which perhaps has a different energy distri¬ 
bution than the 2 + resonance we investigated in the present calculations. We 
find a distribution where the two-neutron virtual .s-state dominates whereas 
the measurements are broader as expected from non-resonance decays. In this 
work we focus on the decay of ” populated” resonances, and an appropriate 
description of this reaction goes beyond the scope of the paper. 

An attempt to understand the distribution was published soon after the ex¬ 
periment in [4], The measured distribution was fitted by a linear combination 
of the lowest hyperharmonic functions of K = 2 and 4. The conclusion was 
that a substantial K = 4 component is needed to reproduce the experiment. 
The decay mechanism dominated by the neutron-neutron virtual s-state was 
abandoned in favor of the K = 4 component. This phenomenological analy¬ 
sis provides a good fit even for energies above the maximum allowed by the 
resonance energy. The decaying resonance wave function does not enter any¬ 
where. The significance is not easy to interpret in terms of decay mechanisms 
as attempted in the present work. 

The neutron energy distribution is not measured but for future comparison we 
show our prediction in Fig. 8. The division into different adiabatic components 
show that the broad total distribution centered around an intermediate energy 
is obtained by adding several qualitatively similar contributions. The different 
mechanisms would all produce most likely energies around half the maximum 
value. To distinguish it is therefore necessary to measure both a-particles and 
neutrons after the decay. 


4-3 Dependence on scaling angle 


It is instructive to investigate the dependence of the distributions on the 
choices of p, 9 and basis size. The p-dependence is already indicated in Fig. 6 
where the distributions are very stable as soon as p is larger than about 50 fm. 
However, this stability does require a sufficiently large basis which at least up 
to 100 fm still can be handled in modest-size computers. It is also clear that a 
finite basis cannot accurately describe the solutions when p increases towards 
cx). Then the angular solutions approach the hyper harmonics but a large basis 
is still required to reproduce the structures at small distances between pairs 
of particles. Eventually this becomes impossible. The many basis functions 
cancel each other at larger distances. 
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At intermediate distances, where the basis is sufficiently large, the resonance 
wave functions are independent of p. For the radial solution this is seen in 
the right part of Fig. 4, where for p larger than about 50 fin, the ratio be¬ 
tween the different radial components is approximately constant. The energy 
distributions are mainly dominated by the absolute squares of these ratios, 
although when different adiabatic components interfere also the real parts of 
these complex ratios may contribute individually. This behavior of the radial 
ratios is responsible for the stable behaviour of the energy distributions for 
sufficiently large values of rho. The constant behavior of the radial ratios also 
implies that the radial wave functions have already reached asymptotics as 
given in eq.(22) for all the channels, and therefore the distributions are inde¬ 
pendent of the scaling angle. However, the latter conclusion is based on an 
assumption of analyticity of the angular solutions as function of 6. When this 
scaling angle is changed corresponding to a rotation across a singularity like a 
two- or three-body resonance the continuity is broken and the solutions change 
as well. 

This is especially clear when we compare two solutions with 6 smaller and 
larger than the angle corresponding to a two-body resonance. For the large 
6 one angular eigenvalue changes character and increases towards infinity as 
p 2 , see [26]. This qualitative change of behavior necessarily causes a change of 
the angular wave functions because the upgoing eigenvalue at large distances 
fully describes the properties of the two-body resonance. These features were 
distributed over several wave functions for the small d-value. 

In between singularities the individual angular solutions are independent of 
both p and 6. This may not be an apparent feature of the numerical solutions 
because the basis has to be sufficiently large for a complete description. As 6 
increases the effective ranges of the two-body interactions also increase and 
the stable region is pushed to larger p-values. This means that the minimum 
basis size has to increase with 6. 


5 Summary and conclusions 


We formulate a method to compute the energy distribution of three parti¬ 
cles emerging after three-body decay of a many-body resonance. The complex 
energy of a resonance corresponds to a pole in the momentum-space wave func¬ 
tion which has an absolute square of the form as Breit-Wigner shape multiplied 
by a smoothly varying function. In coordinate-space this form corresponds to 
a large-distance asymptotic wave function consisting of only outgoing waves. 
We show formally by Fourier transformation that the coordinate-space asymp¬ 
totic angular dependence determines the energy distribution by substituting 
momentum directions for the conjugate coordinate directions. For this the 
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divergent Fourier integral is regularized by the Zeldovich prescription. 


For two-body decay the energy distribution is trivially given by the Breit- 
Wigner distribution of the initial resonance. Energy conservation is taking 
care of everything else. For three-body decay the total energy can be dis¬ 
tributed continuously among the three particles. We show that the resonance 
decay results in distributions obtained from the large-distance angular behav¬ 
ior of the coordinate wave function. The asymptotic behavior can correspond 
to either genuine three-body structures or two-body substructures for example 
corresponding to two-body resonances or configurations favored by substantial 
attraction as for virtual states. Also virtual population of two-body interme¬ 
diate substructures is allowed as an appropriate asymptotic behavior with a 
resulting special energy distribution. The different asymptotics characterize 
the different decay modes used in analyses of experimental data. Different 
modes can co-exist. 


We illustrate by the decay of the 2 + -state in 6 He. The practical computations 
employ the hyperspherical adiabatic expansion combined with the complex 
scaling method. We discuss how a large hyperradius necessarily must be ac¬ 
companied by a large basis. Convergent results may then be obtained with 
less efforts at moderate hyperradii and moderate basis sizes. For convergence 
it is crucial to have all three Faddeev components, and especially if all de¬ 
cay mechanisms simultaneously should be included in the theoretical formu¬ 
lation. The wave function undergoes dramatic changes from small distances, 
where the resonance properties usually are determined, and large distances 
from which the energy distributions emerge. The reasons for these structural 
changes are that the small distance behavior is determined by the two-body 
resonance substructures, whereas the large-distance behavior is determined as 
a competition between two effects, i.e. the energetically favored configuration 
of smallest two-body angular momentum with attractive two-body potentials 
and maintaining the same structure as at small distances in higher-lying levels. 


In conclusion, theoretical interpretation of the simplest nuclear three-body 
decay without Coulomb interactions is already rather complicated. It is then 
advisable to test any given method on these systems. The accuracy of compu¬ 
tations of the more complicated decaying charged systems can then be judged. 
This is important since almost all nuclear three-body resonance decays involve 
charged systems. The goal is to interpret the soon-to-come accurate experi¬ 
mental correlation data for three-body decays of charged systems. 
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